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ABSTRACT 

A simplified framework is presented for assessing the qualitative sensitivities of computed microwave 
properties, satellite brightness temperatures, and radar reflectivities to assumptions concerning the physical 
properties of ice-phase hydrometeors. Properties considered included the shape parameter p of a gamma size 
distribution and the melted-equivalent mass median diameter D 0 , the particle density, the dielectric mixing 
formula, and the choice of complex index of refraction for ice. These properties are examined for selected 
radiometer frequencies of 18.7, 36.5, 89.0, and 150.0 GHz and radar frequencies at 2.8, 13.4, 35.6, and 94.0 
GHz — consistent with existing and planned remote sensing instruments. Passive and active microwave ob- 
servables of ice particles are found to be extremely sensitive to the D 0 of the size distribution. Similar large 
sensitivities are found for variations in the ice volume fraction whenever the geometric mass median diameter 
exceeds approximately Vsth of the wavelength. At 94 GHz the two-way path-integrated attenuation is po- 
tentially large for dense/compact particles. The distribution parameter p has a comparatively weak effect on 
any observable: less than 1-2 K in brightness temperature and a maximum of 2.7 dB (S band only) in the 
effective radar reflectivity. Reversal of the roles of ice and air in the Maxwell Garnett dielectric mixing 
formula leads to a substantial change in both microwave brightness temperature (—10 K) and radar re- 
flectivity (approximately 2 dB across all frequencies). The choice of the complex index of refraction of ice can 
produce a 3%-4% change in the brightness temperature depression. 


1. Introduction 

Whether passive or active, microwave-based precipi- 
tation retrieval algorithms depend on accurate quantita- 
tive relationships between precipitating-cloud properties 
and microwave observables (Panegrossi et al. 1998; Petty 
2001; Kulie et al. 2010). Because of a lack of direct 
knowledge of the physical and associated radiative 
properties of the precipitation particles, however, physi- 
cal models used to derive these relationships have his- 
torically relied on a variety of subjective assumptions 
about the microphysical and microwave radiative transfer 
properties of individual hydrometeors within a cloud, for 
example, the shape, size, and composition of raindrops, 
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snowflakes, aggregates, and graupel particles (Savage 
1978; Smith et al. 1992; Matrosov 1998; Matrosov et al. 
2005; Katsumata et al. 2000; Meneghini and Liao 2000; 
Zhao and Weng 2002; Donovan et al. 2004; Surussavadee 
and Staelin 2006; Casella et al. 2008). Moreover, even the 
three-dimensional distributions of bulk hydrometeor con- 
centrations obtained from state-of-the-art cloud-resolving 
models appear to be alarmingly sensitive to the choice of 
microphysical scheme employed (Hashino and Tripoli 
2007; Han et al. 2010; Wu and Petty 2010). 

Figure 1 summarizes the chain of observations, models, 
and physical assumptions necessary to predict satellite- 
observed microwave brightness temperatures or radar 
reflectivities associated with a particular precipitating-cloud 
system and/or surface precipitation rate. Inappropriate 
assumptions at any point in this chain introduce errors 
into model calculations of the microwave observables 
associated with a given environmental state. Whether 
these errors are tolerable depends on the magnitude and 
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Inputs, Methods, & Assumptions Outputs 



Fig. 1. The relationship between key model assumptions and computed microwave observables associated with 

a modeled precipitating-cloud system. 


spectral dependence of the error and on whether the 
physical model is being used solely as theoretical guidance 
in algorithm design or is being depended upon to supply 
the quantitative basis [e.g., radar reflectivity-rainfall rate 
( Z-R ) relationship or brightness temperature database] 
for the retrievals. 

This paper focuses solely on the sensitivity of computed 
passive and active microwave observables to plausible 
variations in the assumed properties of ice -phase hydro- 
meteors in a precipitating-cloud system. Of potential 
concern are 

(i) errors due to the assumed shapes of the individual 
ice particles (e.g., spheres vs more-realistic crystal 
or aggregate shapes), 

(ii) errors due to the assumed size and/or density 
distributions of the particles within a volume of 
air, and 

(iii) errors due to the dielectric constant assumed for 
pure ice and/or for subwavelength-scale mixtures 
of ice and air. 

This paper is concerned strictly with the first and second 
categories of errors above. We note, however, that the 
microwave properties of realistic nonspherical snow and 
ice particles are the subject of ongoing investigation, as 
exemplified recently by Skofronick- Jackson et al. (2008), 
Noh et al. (2006), Hong (2007a, b), Kim et al. (2007, 2008), 


Ishimoto (2008), Liu (2008), Kulie and Bennartz (2009), 
and Petty and Huang (2010). Nonspherical particles, 
while having apparently more “realistic” shapes, are still 
idealized models, the variations of which have not been 
fully examined in the literature. The authors have already 
made simulations of nonspherical particles in specific 
contexts — see, for example, Petty and Huang (2010) for 
Ku- and Ka-band radar properties. 

Notwithstanding recent progress by the present au- 
thors and others listed above, published simulations of 
nonspherical particles are not extensive and varied 
enough to permit systematic sensitivity studies encom- 
passing reasonable variations in particle size distribu- 
tions (PSDs), particle densities, and other relevant 
properties. For example, Petty and Huang (2010) have 
found that a limited set of nonspherical particles is 
generally unable to reproduce the same frequency- 
dependent scattering and extinction properties as spheres, 
no matter how one adjusts the densities. Of specific 
relevance to this work, however, is that computation 
using variable-density spheres “brackets” the currently 
computed nonspherical particle properties for a wide 
variety of frequencies, shapes, and sizes; in this sense, 
then, the study presented here provides a “worst case” 
scenario with respect to uncertainties arising out of 
variations in the simulated physical properties of 
hydrometeors. 
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Furthermore, there is still considerable present-day 
modeling and algorithm-development work that relies 
on the soft-sphere model for lack of good practical 
alternatives. For example, the Goddard profiling al- 
gorithm (Kummerow et al. 1996) that is used in the 
Tropical Rainfall Measuring Mission, Advanced Mi- 
crowave Scanning Radiometer for Earth Observing 
System (AMSR-E), and Global Precipitation Measure- 
ment (GPM) passive microwave retrieval algorithms 
continues to employ spherical particles in its brightness 
temperature computations. 

To summarize, the “soft sphere” model discussed 
herein must be understood as providing useful guidance 
(as opposed to quantitative predictions) on the potential 
magnitude of uncertainties in forward radiative transfer 
models and retrieval algorithms. In the next section, we 
review the conceptual and mathematical framework for 
computing microwave and bulk physical properties from 
assumed distributions of hydrometeors, including spe- 
cific assumptions about the hydrometeor dielectric 
model, bulk density, and size distribution. We specifically 
confine our attention to spherical ice particles having ei- 
ther a constant or a power-law bulk density and obeying 
a gamma PSD. We also discuss choices bearing on the 
assumed dielectric constant of an ice-air mixture, in- 
cluding both the dielectric constant of pure ice and the 
multiple possible “effective medium” approximations — 
in particular, those of Bruggeman (Bruggeman 1935) and 
the Maxwell Garnett methods (Maxwell Garnett 1904). 
We also consider two values for the complex dielectric 
constant of pure ice: one obtained from the widely used 
tabulation of Warren (1984, hereinafter W84) and the 
other from the recent update by Warren and Brandt 
(2008, hereinafter WB08). We conclude that section with 
a simplified two-stream radiative transfer framework for 
evaluating the approximate sensitivity of observed mi- 
crowave brightness temperatures to variations in the 
hydrometeor model employed. 

In section 3, we evaluate the sensitivities of bulk 
hydrometeor radiative properties, satellite microwave 
brightness temperatures, and radar backscatter at se- 
lected frequencies to plausible variations in the mass 
median diameter D 0 , the size distribution shape pa- 
rameter /i, the particle density p, the choice of dielectric 
constant for pure ice, and the choice of dielectric 
mixing formula. For radiometric calculations, the 
chosen frequencies are 18.7, 36.5, 89.0, and 150 GHz. 
For radar backscatter calculations, we consider 2.8, 
13.4, 35.6, and 94 GHz. In the final section, we con- 
clude with a summary of key results from section 3 and 
implications for the retrieval of precipitation and hy- 
drometeor properties from passive and active micro- 


2. Physical and radiative properties 

Microwave radiances and reflectivities observed from 
either space- or surface-based sensors depend on the three- 
dimensional distribution of local radiative transfer prop- 
erties and temperature along and near the line of sight. 
The frequency-dependent and hydrometeor-dependent 
properties of interest here include the volume or mass 
extinction coefficient ( k e or K e ), single-scattering albedo 
(b, scattering asymmetry parameter g, and radar reflec- 
tivity 77 . These quantities are defined for a volume of air 
containing a suspension of particles; they are derived in 
turn from knowledge of the number and radiative cross 
sections and scattering phase functions of the individual 
particles in the volume. 

In this section, we specify a model for the density of 
each particle as a function of particle diameter or mass. 
The density in turn maps to a microwave-frequency- 
dependent dielectric constant according to either of two 
models for the dielectric constant for pure ice and one of 
three models for the dielectric constant of a mixture of 
ice and air. We calculate radiative cross sections for in- 
dividual particles using standard Mie theory (Mie 1908) 
and then integrate over a specified PSD to obtain bulk 
radiative properties for an ensemble of particles. The 
details of these steps follow below. Because we are pri- 
marily interested in the sensitivity of radiative transfer 
calculations to variations in the physical properties of 
ice-phase hydrometeors, we do not consider the effects 
of absorption or emission by atmospheric gases or cloud 
liquid water, nor do we consider variations of surface 
emissivity. 

a. Particle density 

To assess the impact of ice -particle composition on 
the above-described radiative properties, we assume 
homogeneous dielectric spheres of mass M and geo- 
metric diameter D g that are composed of a blend of ice 
and air — with / ice giving the corresponding volume 
fraction of ice such that / ice + / a i r = 1. The bulk particle 
density is then 

p = 6MI(ttDI) = / lceP , = (DID/p, , (1) 

where we take p t = 917 kg m -3 and neglect the mass of 
air (but not the volume) in these mixtures. Because of 
the multitude of morphological possibilities, the re- 
lationship between the geometric size of a snowflake 
and its effective density is not well characterized, but 
it is conventional to model the relationship as a power 
law: 

p( D g ) = a p D g bp . 


wave sensors. 


( 2 ) 
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Table 1. Snow density parameters (centimeter-gram-second 
units) for Eq. (2) from Brown and Ruf (2007). 


Source 

a p 

b P 

Magono and Nakamura (1965) 

0.022 

1.5 

Barthazy et al. (1998) 

0.018 

0.8 

University of Wisconsin Nonhydrostatic 
Modeling System (Tripoli 1992) 

0.015 

0.6 

Locatelli and Hobbs (1974) 

0.015 

1.18 

Mitchell et al. (1990) 

0.012 

1.0 


We examine two cases that commonly occur in the 
literature: 1) the assumption of a constant density 
throughout the PSD ( b p = 0) and 2) a power-law density- 
size relationship, of which five published examples are 
taken from Brown and Ruf (2007) and Szyrmer and 
Zawadzki (2010) and summarized in Table 1. 

b. Effective dielectric constant 

In the spherical-particle approximation, low-density 
ice particles are modeled as having a uniform dielectric 
constant appropriate to a particular mixture of ice and 
air. Even once this mixture has been specified, how- 
ever, a surprising number of choices exist for comput- 
ing the effective dielectric constant. One objective of 
this study is therefore to elucidate the practical sensi- 
tivity of microwave observables to the chosen mixing 
rule. 

The two best-known mixing rules are the Bruggeman 
(Bruggeman 1935) or effective-medium approximation 
and the Maxwell Garnett (Maxwell Garnett 1904) or 
matrix-inclusion method. Both methods continue to be 
widely used in physical modeling of the microwave ra- 
diative properties of hydrometeors. These methods are 
summarized in Bohren and Huffman (1983). Casella 
et al. (2008), Petty and Huang (2010), and others have 
examined the consistency of these two common mixing 
rules with explicit discrete dipole approximation (DD A) 
calculations of mass-weighted extinction, scattering, and 
backscattering properties for idealized reduced-density 
spherical ice particles. 

Among other things, Petty and Huang (2010) found 
that a generalized mixing rule falling roughly midway 
between the two cases gave the best fit to the DDA- 
calculated properties. In this study we extend their 
approach to incorporate a wider range of microwave 
frequencies and particle densities and to examine the 
sensitivity of variations in the previously described 
physical properties of hydrometeors on simulated/ 
idealized passive microwave brightness temperatures 
and radar reflectivities. 

The Bruggeman formulation for a mixture of two 
components is 


e eff e B 2 ^ 2 + 2e l e 2) 2^’ ^) 

where 

jS = l[(l-3/ 2 )£ 2 + (l-3/ 1 )e 1 ]. (4) 


In Eq. (3), e 1 and e 2 are the dielectric constants of 
components 1 and 2, respectively; fi and f 2 = 1 — f\ are 
the respective volume fractions; and e B is the effective 
dielectric constant of the mixture. 

The Maxwell Garnett (MG) method was derived for 
a two-component mixture consisting of a matrix with 
dielectric constant e m and a fraction / of small spherical 
inclusions with dielectric constant e t . The effective di- 
electric constant according to this method is given by 


e eff e MG 


e ( .(l + 2/) — e m (2/ — 2) 

m e m (2+f) + e i (l-f) ' 


(5) 


Reversing the roles of e u e m , and / results in a different 
value for e e ff i n the MG case. 

Figure 2 illustrates the influence of ice-air volume 
fraction on the real and imaginary components of the 
dielectric function for the Bruggeman and two MG 
formulations. The top and bottom panels of Fig. 2 show 
the dielectric permittivities of ice from W84 and WB08, 
respectively; both are presented at 89 GHz and a tem- 
perature of 266 K. Each circle is labeled according to the 
volume fraction of ice. 

It is apparent that at small volume fractions of ice the 
Bruggeman method is closer to MG- ([I] A) , in which ice is 
the inclusion component. As the volume fraction of ice 
increases, however, the Bruggeman curve approaches the 
MG-([A]I) curve, in which air is the inclusion component. 
The two MG methods effectively bracket the Bruggeman 
results (for all frequencies and ice-volume fractions 
tested; not shown), with Bruggeman being closest to the 
MG method where the volume fraction of the inclusion 
components is small (e.g., / inc i < 0.3; Sihvola and Kong 
1988, Lakhtakia and Shanker 1993, Mackay 2005). For 
particles having large size parameters, the size parame- 
ters of the individual inclusions may become unphysically 
large for the MG methods to be appropriate; therefore, 
additional care must be taken to ensure that this does not 
occur if one uses the MG method. 


c. Single-particle properties 

Any single hydrometeor can be characterized physi- 
cally by its mass M, size, shape, and composition (e.g., 
fraction of ice vs liquid). Of these properties, only the 
mass and composition are unambiguous. Geometric 
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Permittivity of ice-air mixture at 89 GHz, T = 266K 
[Warren, 1984] 



Permittivity of ice-air mixture at 89 GHz, T = 266K 
[Warren and Brandt, 2008] 



Fig. 2. Dielectric permittivity of ice-air material mixture at 
89 GHz and 266 K for ice volume fractions ranging from 0.01 to 
1.0. The curves represent the Maxwell Garnett and Bruggeman 
formulations. Label MG- ([A] I) indicates an ice matrix with air 
inclusions; MG- ([I] A) denotes the reverse, (top) The W84 per- 
mittivity; (bottom) the permittivity from WB08. 




/6 M \ 1/3 

WpJ 


where p w is the density of water, taken to be 
1000 kg m“ 3 . 

For a given microwave frequency, each particle also 
has an absorption cross section cr a , scattering cross sec- 
tion cr s , and scattering phase function P(0). For true 
nonspherical particles, these characteristics would ex- 
hibit a strong dependence on orientation; such effects 
fall outside the scope of this paper. 

Given the geometric diameter D g and complex index 
of refraction N, where N 2 = e e ff> the assumption of 
spherical shape allows us to use standard Mie theory 
(Mie 1908; Bohren and Huffman 1983) to compute the 
above properties and, from these, the extinction cross 
section: 


and the radar backscatter cross section: 

^=‘U > (©)I @=17 - (7) 

d. Bulk microphysical properties 

Atmospheric remote sensing methods respond not to 
individual particles but rather to large numbers of sus- 
pended particles of different sizes, shapes, and orienta- 
tions distributed more or less randomly throughout 
a substantial volume of the atmosphere. When dealing 
with single-phase ice particles, it is convenient (although 
by no means strictly accurate) to assume that the sole 
relevant independent variable is the particle mass M or, 
alternatively, the melted-equivalent diameter D. 

Size distribution 

Following Ulbrich (1983), Brandes et al. (2004), and 
Heymsfield et al. (2005), and many others, we assume 
a three-parameter gamma PSD relative to the melted- 
equivalent diameter D : 

n(D ) = N 0 D^ exp(- A D ) , (8) 

where A is often referred to as the “slope parameter,” p 
is a “shape parameter,” and N 0 controls the overall 
concentration of particles. 1 * It is important to note that 


size, shape, and effective density are far more variable 
and are difficult even to define, especially in the case of 
complex snow aggregates (Petty and Huang 2010, 2011). 
One may also characterize the mass of an ice particle by 
its melted-equivalent spherical radius D, defined as 


1 As discussed by Petty and Huang (2011), a gamma distribution 
in D is also a gamma distribution in D g only when density is con- 

stant. For the case in which the density is nonconstant and given by 

Eq. (2), the corresponding PSD expressed in D g is a four-parameter 
modified gamma distribution. 
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the use of a continuous gamma distribution is not ideal 
for data for which the tails of the distribution may be 
underrepresented, such as in observational datasets. In 
such cases we recommend using a truncated gamma 
distribution (Ulbrich and Atlas 1998; Hogan et al. 2012). 
Here, we prefer the continuous gamma distribution for 
computational convenience. 

The mass median melted-equivalent diameter D 0 is 
the characteristic size of the PSD that divides the total 
mass of the PSD into two equal parts and is given for the 
gamma distribution by 


cos0) of the ensemble scattering phase function P(0) 
(Bohren and Huffman 1983). The latter is in turn the 
scattering cross-section-weighted sum of the phase 
functions of the individual particles. 

In the absence of attenuation, the effective radar 
reflectivity factor Z e is given by 

A 4 f 00 

a b (D)n(D)dD, (13) 

^ \ K J 

where A is the radar wavelength and \K W \ 2 = 0.93 by 
convention. 

f Radiative transfer 


The corresponding mass median geometric diameter 
D 0g is obtained with Eq. (1), and the ice water content is 
given by 



( 10 ) 


e. Bulk radiative properties 

The radiative cross sections of individual particles are 
computed for the microwave frequencies of interest by 
using an updated version of the Mie routine published 
by Bohren and Huffman (1983). They are then summed 
over the PSD to yield the total cross sections per unit 
volume of air and related radiative properties. Of par- 
ticular relevance are the volume extinction coefficient, 


k = 


*00 

o- (D)n(D)dD, 

Jo 


( 11 ) 


and the volume scattering coefficient k s , defined simi- 
larly using the scattering cross section cr s {D). The en- 
semble single-scatter albedo is co = k s /k e . 

Here it is convenient to characterize extinction not in 
terms of the volume extinction k e (m -1 ) but rather in 
terms of the mass extinction coefficient K e (m 2 kg -1 ): 

= ( 12 ) 


which represents the extinction cross section per unit 
mass along the line of sight. It has two advantages in the 
present context: first, it depends only on the shape of 
n(D ) and not on the magnitude; second, it allows optical 
path and transmission to be conveniently related to mass 
path along the line of sight rather than to the geometric 
distance. 

For most radiative transfer calculations, we addi- 
tionally require at least the asymmetry parameter g, 
which is defined as the first moment (with respect to 


The purpose of this paper is not to undertake accu- 
rate and detailed radiative transfer calculations but 
rather to evaluate relative sensitivities of microwave 
observables to varying assumptions about ice particle 
properties their PSDs. We may therefore take advan- 
tage of highly simplified treatments of the radiative 
transfer to characterize the expected magnitude of these 
responses. 

1) Brightness temperature 

The general theoretical basis for computing radiative 
transfer in a scattering and emitting medium is well 
known and need not be repeated here (Lenoble 1985; 
Liou 2002). For the purpose of sensitivity analysis, we 
introduce a simple analytic approximation that is based 
on a two-stream approximation (TSA) that has been 
found to be remarkably accurate at describing thermal 
emission from a plane -parallel slab viewed at the —55° 
incidence typical of a satellite microwave imager. Weng 
and Grody (2000) and Liu and Weng (2002) have also 
discussed the utility of the TSA in a similar context; 
here, we develop expressions that starting with the TSA 
formulation as reviewed by Petty (2006). 

For an isothermal finite slab of hydrometeors with 
temperature T overlying a nonreflective surface with 
temperature T s , the total optical thickness is t = K e l , 
where X is the ice water path, and the upwelling mi- 
crowave brightness temperature in the two-stream ap- 
proximation is 


T B = T s t+(l-r-t)T, (14) 

where t is the total transmittance (direct plus diffuse) 
and r is the reflectance. In the two-stream approxima- 
tion, the reflectance and transmittance are in turn 
functions of r and of the bulk single-scattering albedo co 
and asymmetry parameter g of the constituent particles 
within the slab: 
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where 


rje_2l lT ) 

p Yt — r 2 p~Tr 


and 


t = - 


1 -r 2 

1 'oc 


_ A A -g 

^/l -g&»+ V / U r g 

r^2[(i-«g)(i-(w)] 1/2 


and 


(15) 

(16) 


(17) 

(18) 


2) Effective radar reflectivity 

At the top of the slab, Z e is calculated using Eq. (13) 
applied to the assumed particle distribution. At the 
bottom of the slab, Z e is reduced by the two-way path 
attenuation. We ignore multiple scattering in these 
idealized sensitivity experiments, although we note that 
Matrosov et al. (2008) and Battaglia et al. (2008) have 
shown that multiple scattering (MS) can significantly 
enhance the observed reflectivity profiles at 94 GHz 
above what single scattering alone would predict. At 
13.4 and 35.6 GHz, the reflectivities appear to be rela- 
tively unaffected by MS in snowfall observations. 

For nadir viewing by a satellite radar, the direct trans- 
mittance to the base of the slab is given by Beer’s law: T = 
exp(-r) and may be expressed as a one-way path atten- 
uation in decibels with A (dB) = -10 logi 0 T. The two- 
way attenuation is, of course, 2 times this amount. 

We may thus not only examine the effect of assumed 
particle properties on the intrinsic reflectivity Z e but also 
on the combined effect of changes in reflectivity and 
changes in two-way attenuation: 

z; (dBZ) - z e (dBZ) - 2A (dB) . (19) 

Unlike the case for radiometric quantities, which de- 
pend only on the path-integrated hydrometeor properties, 
reflectivity also depends on hydrometeor concentration 
per unit volume. Here, we assume a mass concentration 
W of 1 g m -3 . For a typical mass-weighted fall speed of 
lms -1 , this value would imply a liquid-equivalent pre- 
cipitation rate of 3.6 mm h -1 . 


3. Sensitivity studies 

a. Overview 

The variable parameters of p [Eq. (1)], /x [Eq. (8)], and 
D 0 [Eq. (9)] fully characterize the composition and 
distribution of the ice-phase hydrometeor ensemble 


apart from the multiplicative constant N 0 . In this paper, 
N 0 is required only for computing radar reflectivity Z e . 
We set Nq implicitly by specifying a fixed ice water 
content W of 1 g m -3 , unless otherwise noted. 

Choosing the microwave frequency v and one of three 
candidate dielectric mixing formulas then allows all 
local microwave properties K e , co, g, and Z e to be com- 
puted, along with T B and A following from the re- 
lationships in section 2f and from the specified column 
ice water path X of 1 kg m -2 . 

We select four frequencies — 18.7, 36.5, 89.0, and 150.0 
GHz — that are representative of those found on current 
and planned microwave imagers. At each frequency, we 
examine the effects of varying p, D 0 , p, and the mixing 
formula on the computed single-layer slab TSA T B 
values. We also performed calculations for 183.3 GHz 
and found that it behaves qualitatively very much like 
150 GHz; general conclusions in this paper concerning 
the latter frequency also apply to 183.3 GHz. 

For our examination of the radar reflectivities Z e and 
A , we select an additional four frequencies: 2.8, 13.4, 
36.5, and 94.0 GHz. The first of these corresponds to the 
frequency utilized by the ground-based Weather Sur- 
veillance Radar-1988 Doppler (WSR-88D). The next 
two will be utilized by the GPM dual-frequency pre- 
cipitation radar (DPR). The 94.0-GHz frequency is 
utilized by the CloudSat radar. 

For each of the above frequencies, we computed the 
relevant radiative properties by assuming spherical 
particles having either a constant density ( b p = 0.0) or 
a power-law density defined by Eq. (2) and Table 1. We 
vary p (through a p and b p ), D 0 , and p for each of the 
three dielectric mixing formulas and for both the WB08 
and older W84 tabulations of the complex index of re- 
fraction of ice. 

Representative values of all of these parameters 
were chosen to define a baseline against which all var- 
iations were compared. For this baseline, we assume a 
constant-density PSD with / ice = 0.1 (p ~ 91.7 kg m -3 ), 
D 0 = 1 mm, and p = 0 [i.e., n(D) is exponential]; e ice is 
given by WB08, and the dielectric mixing rule is 
Bruggeman (B). 

One parameter at a time was varied from the baseline 
value. In this case, D 0 was varied from 0.125 mm, cor- 
responding to fairly small cirrus particles, to 8.0 mm, 
corresponding to large sleet, snow aggregates, or small 
hail. The ice fraction / ice was varied from 0.01, corre- 
sponding to a very-low-density snow aggregate, to 1.0, 
corresponding to a solid ice pellet. The shape parameter 
p was varied from - 1 to 4, encompassing the majority of 
empirically determined values for this parameter (e.g., 
Heymsfield 2003). Representative results are tabulated 
in Tables 2-4, described in more detail below. 
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Table 2. Local radiative transfer properties of varying ice particle distribution models for selected radiometer frequencies. Baseline 
parameter values are indicated in boldface. For ease of comparison, baseline results are repeated within each parameter grouping. The 


square brackets under “mixing rule” enclose the Maxwell Garnett inclusion component, either ice or air. Note that 

K e has units of meters 

squared per kilogram. 










36.5 GHz 



89.0 GHz 



150 GHz 



«e 

QJ 

g 

K e 

O) 

8 

K e 

O) 

8 

/ice (constant p) 

0.01 (0.009) 

0.80 X 10“ 2 

0.896 

0.68 

0.59 X 10" 

_1 0.917 

0.90 

0.18 

0.921 

0.96 

0.03 (0.028) 

0.13 X 10" 1 

0.937 

0.52 

0.11 

0.955 

0.84 

0.35 

0.958 

0.92 

0.1 ( 0 . 092 ) 

0.23 X 10" 1 

0.960 

0.33 

0.24 

0.976 

0.73 

0.78 

0.978 

0.86 

0.3 (0.28) 

0.38 X 10' 1 

0.971 

0.21 

0.51 

0.985 

0.59 

1.7 

0.986 

0.76 

1.0 (0.92) 

0.48 X 10" 1 

0.980 

0.15 

1.00 

0.990 

0.44 

3.1 

0.986 

0.53 

a p , b p [Eq. (2)] 
0.022, 1.50 

0.75 X 10" 1 

0.964 

0.26 

1.12 

0.983 

0.50 

4.86 

0.985 

0.65 

0.012, 1.00 

0.26 X 10' 1 

0.961 

0.32 

0.27 

0.977 

0.72 

0.83 

0.979 

0.85 

0.015, 1.18 

0.31 X 10" 1 

0.960 

0.32 

0.34 

0.977 

0.68 

1.16 

0.980 

0.80 

0.018, 0.80 

0.25 X 10' 1 

0.961 

0.32 

0.26 

0.977 

0.72 

0.87 

0.979 

0.85 

0.012, 0.60 

0.24 X 10" 1 

0.961 

0.32 

0.25 

0.977 

0.72 

0.83 

0.979 

0.86 

D 0 (mm) 
0.125 

0.99 X 10“ 3 

0.091 

0.01 

0.83 X 10" 

" 2 0.356 

0.04 

0.36 X 10 _1 

0.571 

0.12 

0.25 

0.16 X 10' 2 

0.433 

0.03 

0.24 X 10" 

_1 0.774 

0.16 

0.11 

0.860 

0.35 

0.5 

0.56 X 10“ 2 

0.836 

0.11 

0.85 X 10" 

_1 0.934 

0.43 

0.33 

0.950 

0.66 

1.0 

0.23 X 10' 1 

0.960 

0.33 

0.24 

0.976 

0.73 

0.78 

0.978 

0.86 

2.0 

0.73 X 10" 1 

0.987 

0.65 

0.56 

0.989 

0.89 

1.58 

0.989 

0.95 

4.0 

0.18 

0.995 

0.86 

1.08 

0.994 

0.96 

2.48 

0.993 

0.97 

8.0 

0.37 

0.997 

0.94 

1.54 

0.996 

0.97 

2.39 

0.993 

0.89 




-1.0 

0.26 X 10" 1 

0.963 

0.39 

0.25 

0.977 

0.75 

0.81 

0.979 

0.87 

0.0 

0.23 X 10" 1 

0.960 

0.33 

0.24 

0.976 

0.73 

0.78 

0.978 

0.86 

1.0 

0.22 X 10" 1 

0.957 

0.30 

0.23 

0.975 

0.71 

0.77 

0.978 

0.86 

2.0 

0.21 X 10" 1 

0.956 

0.27 

0.23 

0.975 

0.70 

0.76 

0.978 

0.85 

3.0 

0.20 X 10" 1 

0.954 

0.25 

0.22 

0.974 

0.69 

0.75 

0.978 

0.85 

4.0 

0.20 X 10" 1 

0.953 

0.23 

0.22 

0.974 

0.69 

0.75 

0.977 

0.85 

: ice 

WB08 

0.23 X 10" 1 

0.960 

0.33 

0.24 

0.976 

0.73 

0.78 

0.978 

0.86 

W84 

Mixing rule 

0.25 X 10" 1 

0.894 

0.33 

0.24 

0.953 

0.73 

0.78 

0.966 

0.86 

MG-([I]A) 

0.22 X 10" 1 

0.961 

0.33 

0.22 

0.977 

0.73 

0.73 

0.979 

0.86 

Bruggeman 

0.23 X 10" 1 

0.960 

0.33 

0.24 

0.976 

0.73 

0.78 

0.978 

0.86 

MG-([A]I) 

0.37 X 10" 1 

0.955 

0.34 

0.38 

0.973 

0.73 

1.2 

0.975 

0.86 


b. Effects on local optical properties 

A number of generalizations may be made about the 
effects of model parameters on the fundamental micro- 
wave radiative properties K e , m, and g shown in Table 2. 
Note that the 18.7-GHz results discussed below are not 
shown in Table 2. 

First, it is apparent that, for fixed D 0 = 1 mm, the as- 
sumed particle density has a profound effect on the mass 
extinction coefficient K e at all frequencies. The higher the 
frequency is, the more dramatic is the effect. In par- 
ticular, the difference in mass extinction coefficient 
between/i ce = 0.01 and / ice = 1 (constant-density PSDs) 
ranges from slightly more than a factor of 2 at 18.7 GHz 
to a factor of 17 at 150 GHz. 

For the variable-density relationships, when b p is 
greater than unity, there is a significantly stronger 


sensitivity in K e than for when b p is less than unity; over 
the five examples, we find as much as a factor-of-6 
difference at 150 GHz. At all frequencies, the setting of 
a p = 0.022 and b p = 1.50 produces larger K e values than 
do even solid ice spheres. 

Higher-density spheres have higher single-scatter al- 
bedos and smaller values of g. The latter effect is ex- 
pected because high density implies a smaller geometric 
dimension relative to the wavelength and thus more 
symmetric scattering. In particular, in the Rayleigh limit 
(A* <C A), we expect g = 0. At the other extreme, low- 
density spheres exhibit very strong forward scattering, 
with g approaching the theoretical limit of 1 at the 
highest frequency. 

Not surprising is that variations in D 0 also strongly 
affect microwave properties, with larger particle masses 
being associated with sharply higher extinction per unit 
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mass. The sensitivity is largest at the lowest frequencies. 
A factor-of-4 change in D 0 (from 0.5 to 2.0 mm) gives 
rise to a factor-of-19 difference in K e at 18.7 GHz but 
only a factor-of-5 difference at 150 GHz. Again, this 
spectral dependence is not unexpected, because the 
transition from lower to higher frequency is associ- 
ated with a move farther away from the Rayleigh re- 
gime, with its D 6 scattering dependence, into the Mie 
regime, where extinction efficiencies are more nearly 
constant. 

The distribution shape parameter p, turns out to have 
only a very subtle effect on radiative properties. For 
fixed D 0 and fixed W, increasing values of p act to 
sharpen and narrow the particle mass distribution 
D 3 n(D ) — that is, there is an increase in the mass con- 
tribution of particles with sizes near D 0 with compen- 
sating decreases at the large and small ends of the 
distribution. At the lowest frequency, K e is approximately 
50% larger for p = -1 than it is for p = 4. 

The choice of W84 or WB08 to specify the complex 
index refraction of ice primarily affects the imaginary 
part, which is roughly a factor-of-2 smaller in the newer 
compilation. This difference gives rise to a roughly 
proportional change in the absorption coefficient of the 
particles and thus in the single-scatter coalbedo 1 — a). 
Where the latter is moderately large (primarily at the 
lowest frequencies), there is also a noticeable increase in 
the overall extinction coefficient. 

Furthermore, we consider the choice of the dielectric 
mixing formula. The baseline choice is the Bruggeman 
formulation (B). The two alternatives considered are the 
Maxwell Garnett formulation with ice as the inclusion 
(MG-([I]A)) or with ice as the matrix (MG-([A]I)). We 
see that B and MG- ([I] A) yield very similar results 
overall while markedly higher absorption and overall 
extinction are associated with MG- ([A] I), particularly at 
the lowest frequencies. This effect has been noted in 
previous studies, and we corroborate it here as well. 

In view of the very strong dependence of K e on both 
D 0 and p, it is worth considering the full range of com- 
binations of these two variables as in Fig. 3. These 
contour plots reveal that the maximum in K e is achieved 
for very dense particles with D 0 ^ A. From this maxi- 
mum, K e falls off by at least two orders of magnitude for 
D 0 < 0.5 mm. As frequency increases, so does the sen- 
sitivity to particle density, which we suspect will be true 
with nonspherical particles as well, providing a basis for 
determining particle composition/density if the charac- 
teristic sizes of the particles are known (e.g., from dual- 
wavelength radar observations). 

In Fig. 4a, we have plotted K e versus D 0 for the five 
selected size-density relationships (points). Also in- 
cluded, for reference, is K e for selected constant-density 


values (curves). It is evident that the variable-density 
quantities exhibit striking differences relative to the 
constant-density curves, especially in the two cases in 
which b p is greater than unity. This point is especially 
relevant to researchers who employ constant-density 
spheres in that they may be unnecessarily biasing their 
results toward lower extinction values for a given D 0 (or 
ice water content) value. 

1) Results for unit ice water path 

We have previously stated that the two-stream ap- 
proximation, while novel, is sufficient to qualitatively 
explore the brightness temperature response to varia- 
tions of the physical properties of ice hydrometeors 
present in the atmosphere. With this context in mind, we 
now consider the combined effect of K e , co, and g on 
characteristic brightness temperature depressions rela- 
tive to the 273-K background, using the aforementioned 
X = 1 kg m -2 as the benchmark. For other values of Z, 
brightness temperatures can readily be computed from 
Table 2 and Eq. (14); in general, however, one can ex- 
pect the brightness temperature depression to be 
roughly proportional to X at least down to several tens of 
kelvin below 273 K. 

In addition to the tabulated values in Table 3, we 
examine the response to joint variations in D 0 and p in 
Figs. 4-6. In considering these together, we note the 
following behavior: 

1) For fixed ice water path, the brightness temperature 
depression is very sensitive to both D 0 and p. The 
minimum T B (maximum scattering-induced depres- 
sion) occurs for the variable-density formulation of 
Magono and Nakamura (1965) in Fig. 4c (only 
150 GHz is shown). For constant-density PSDs, 
D 0 ~ 0.6A (Fig. 5), and the strongest T B depressions 
occur at the highest-density particles (solid spheres). 

2) As expected from the results in Table 2, the 
parameter p has a nearly negligible effect on T B 
for modest ice water path, amounting to a maximum 
spread of 1.6 K at 150 GHz when D 0 = 1 mm and 
p = 0.092 g cm -3 (/ ice = 0.1). We examined the 
response for the full range of p and D 0 (not shown) 
and found that p only begins to play a significant 
role in T B at 150 GHz for D 0 > 4 mm, in which case 
T b may be as much as 10-20 K greater for p = 4 
than it is for p = — 1. 

3) The consequences of choosing W84 or WB08 for 
the index of refraction are potentially substantial, 
especially when considering the selected mass-size 
relationships. While the maximum response for base- 
line (p = 0.092 g cm -3 and D 0 = 1 mm) is only a few 
tenths of a degree for any frequency, the full range 
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Mass Extinction Coefficient [m 2 /kg] 

a) 18.7 GHz b) 36.5 GHz 



c) 89.0 GHz d) 150.0 GHz 



Fig. 3. The mass extinction coefficient K e of ice obeying an exponential size distribution depicted as a function of the particle 
density p and melted equivalent mass median diameter D 0 (contour values are logarithmically spaced) for (a) 18.7, (b) 36.5, (c) 89.0, 
and (d) 150.0 GHz. 


of p and Do has a maximum sensitivity of as much as 
14° for 89 GHz and 27° at 150 GHz— mainly for 
denser particles (see Fig. 6). 

4) The choice of dielectric mixing formulations has 
a minor impact at lower frequencies. However, T B 
depressions at 89 and 150 GHz are “artificially” 
enhanced when using MG- ([A] I) because of the 
asymmetric treatment of ice as the matrix material. 


2) Two-channel ratio of T b depressions 

For any given microwave frequency, the scattering- 
induced brightness temperature depression S , defined in 
the context of this analysis as 

S = 273-T B , (20) 
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Spherical Particle Properties @ 150 GHz 
constant density (solid lines) & 5 size-density relationships (points) 





D 0 [mm] D Q [mm] 

Fig. 4. Distribution-averaged properties at 150 GHz: (a) mass extinction coefficient, (b) single-scattering albedo, (c) the two-stream 
brightness temperature associated with an isothermal (X = 266 K) layer of exponentially distributed ice particles (X = 1 kg m -2 ) and 
melted equivalent diameter D 0 overlying a black surface having temperature 273 K, and (d) the scattering asymmetry parameter. Solid 
lines are associated with constant density throughout the PSD; points indicate variable size-density relationships [Eq. (2)]. 


is clearly a sensitive function of three variables: X, Z) 0 , and 
p. It can be shown that, for fixed frequency, particle size, 
and density, S is approximately proportional to X as long as 
the ice-layer optical thickness is less than about unity. But 


since D () and p are unknown and strongly affect S as well, 
single-channel retrievals of X are unlikely to be meaningful. 

By using multiple radiometer frequencies, one might 
hope to isolate the effects of X, D 0 , and p, allowing 
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Table 3. Radiometric response to ice particle distribution pa- 
rameters as shown by the brightness temperature T B (K) for an ice 
water path T of 1 kg m -2 . Ice-layer temperature is 266 K; back- 
ground brightness temperature is 273 K. Baseline parameter 
values are indicated in boldface. For ease of comparison, baseline 
results are repeated within each parameter grouping. 



18.7 GHz 

36.5 GHz 

89.0 GHz 

150 GHz 

/ ice (constant p) 

0.01 (0.009) 

272.8 

272.4 

271.5 

271.1 

0.03 (0.028) 

272.6 

271.4 

268.5 

265.9 

0.1 (0.092) 

272.4 

269.0 

256.9 

247.3 

0.3 (0.28) 

272.2 

265.3 

226.9 

197.5 

1.0 (0.92) 

272.2 

262.5 

177.2 

122.0 

a p , b p [Eq. (2)] 

0.022, 1.50 

271.2 

259.3 

178.5 

117.2 

0.012, 1.00 

272.3 

268.5 

254.2 

242.6 

0.015, 1.18 

272.2 

267.6 

247.5 

224.8 

0.018, 0.80 

272.3 

268.6 

254.5 

243.4 

0.015, 0.60 

272.4 

268.7 

255.3 

245.1 

D 0 (mm) 

0.125 

273.0 

273.0 

272.2 

268.1 

0.25 

273.0 

272.8 

268.8 

257.0 

0.5 

272.9 

271.9 

261.2 

247.3 

1.0 

272.4 

269.0 

256.9 

247.3 

2.0 

270.9 

266.3 

257.4 

253.4 

4.0 

269.5 

266.3 

261.8 

254.6 

8.0 

269.4 

267.1 

261.0 

218.0 

p 

-1.0 

272.3 

268.9 

257.3 

248.2 

0.0 

272.4 

269.0 

256.9 

247.2 

1.0 

272.4 

269.0 

256.2 

247.6 

2.0 

272.5 

269.0 

256.0 

246.3 

3.0 

272.5 

269.1 

255.7 

246.5 

4.0 

272.5 

269.1 

255.9 

246.6 

^ice 

WB08 

272.4 

269.0 

256.9 

247.3 

W84 

272.4 

269.0 

257.1 

247.7 

Mixing rule 

MG-([I]A) 

272.4 

269.2 

258.1 

248.7 

Bruggeman 

272.4 

269.0 

256.9 

247.3 

MG-([A]I) 

272.0 

266.8 

248.6 

236.2 


reasonably unambiguous retrievals of each (Weng and 
Grody 2000; Zhao and Weng 2002). The success of this 
approach depends on there being a significant spectral 
dependence of S on D 0 and/or p. One measure of such 
a dependence is the ratio of scattering depressions S\/S 2 
for pairs of frequencies v 1 and v 2 . 

Figures 7a and 7b depict contours of the depression 
ratios S 89 /S 36 while Figs. 7c and 7d show the contours for 
SisoAW These results indicate that the depression ratio 
is a sensitive indicator of D 0 . In particular, for all but the 
lowest-density particles, Fig. 7a shows that measure- 
ments of S 89 /S 36 allow D 0 to be estimated fairly un- 
ambiguously between approximately 0.1 and 3 mm 
while Fig. 7c shows that Si 5 q/S 89 resolves variations of D 0 
in the range of 0.2-1. 0 mm. 


The sensitivity to particle density p is very low for 
ratios of either pair of frequencies. Yet we have seen in 
Fig. 5 that, for fixed D 0 , the magnitude of the response of 
any of these three frequencies depends strongly on both 
X and p. The variable-density relationships (Figs. 7b and 
7d) illustrate a more complex story. It is evident that 
given a ratio of S 89 /S 36 one can find a range of sensitivity 
to particle sizes between 0.1 and 7 mm, depending on 
which density model is assumed. Comparing this with 
^ 150^89 allows for the selection of which of the two 
possible D 0 values is appropriate (i.e., with which side of 
the peak on is concerned), as long as one sticks with 
a given size-density relationship (or constant density 
value). The implication is that, although retrievals of D 0 
may be straightforward, the separation of X from p may 
not be. Direct retrievals of X will most likely depend on 
reasonably confident specifications of the effective ice 
particle density p. 

c. Effects on observed radar backscatter 

We now turn our attention to the effects of ice-model 
parameters on backscattered microwave radiation as ob- 
served by ground-based or space-based radars. There are 
two major components to this response: the intrinsic ef- 
fective radar reflectivity Z e of the distribution and the 
modification of this apparent reflectivity due to attenua- 
tion by the intervening medium. For a given nonzero mass 
path, it is possible for variations in model parameters to 
alter both Z e and the two-way path attenuation ( '2A ) in 
ways that either reinforce or partially cancel each other. 

For our analysis, we choose four representative fre- 
quencies corresponding to the ground-based WSR-88Ds 
(2.8 GHz), the two channels of the spaceborne GPM 
DPR (13.4 and 36.5 GHz), and the CloudSat profiling 
radar (94 GHz). Sample results are presented in Table 4. 
As before, we also examine the full range of combina- 
tions of D 0 and p for selected frequencies in Figs. 8 and 9. 

Seven general results were obtained: 

1) For 13.4, 35.6, and 94 GHz, the effective reflectivity 
per unit mass is dominated by Rayleigh scattering 
and is nearly independent of p for D 0g (and not D 0 ) of 
less than approximately 3, 2, and 0.5 mm, respec- 
tively. For 2.8 GHz (not shown), the corresponding 
threshold is approximately 14 mm. Within this re- 
gime, single-frequency estimates of ice or snow water 
content require only accurate knowledge of the 
melted-equivalent mass median diameter D 0 , assum- 
ing an ideally sensitive radar. 

2) For D 0g that is much larger than the above thresh- 
olds, Z e depends equally strongly on both D 0 and p. 

3) For 13.4 and 36.5 GHz, attenuation per unit mass is 
much less than 1 dB except for large D 0 combined 


D 0 [mm] 
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TB [K] for 1 kg/m 2 IWP over 273 K background 

a) 18.7 GHz b) 36.5 GHz 





Fig. 5. The brightness temperature associated with an isothermal (T = 266 K) layer of exponentially distributed ice particles (T = 
1 kg m -2 ) with density p and melted equivalent diameter D 0 overlying a black surface having temperature 273 K: (a) 18.7, (b) 36.5, (c) 89, 
and (d) 150 GHz. Note the use of a linear axis for D 0 in (a) and (b) and a logarithmic axis in (c) and (d). 


with moderate to high density, in which case it can be 
as much as 6 dB (kg m -2 ) -1 at 35.6 GHz. 

4) For D () = 1 mm and p = 0.092 g cm -3 (/ ice = 0.1), the 
parameter p, has no more than a few tenths of 
a reflectivity decibel effect on either Z e or A at any 
frequency. Only at large sizes (not shown) does the 
effect of jul start to be significant, and even then it is 
minor relative to the effect of variations in D 0 and p. 


5) The choice of WB08 or W84 for the dielectric 
constant makes no significant difference for radar 
reflectivity or attenuation. 

6) As before, the use of the MG-([A]I) mixing 
model enhances scattering by a couple of decibels 
relative to Bruggeman (B) or MG-([I]A). At 
94 GHz, increased attenuation over 1 kg m -2 is 
also significant. 


D 0 [mm] D 0 [mm] 
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a) 89.0 GHz 


T b (WB08) - T b (W84) [K] 

b) 89.0 GHz 



c) 150.0 GHz 




Fig. 6. Similar to Fig. 5, but contours depict the change in brightness temperature that results from choosing the W84 or WB08 values for 
the dielectric constant of ice: (a) 89 and (c) 150 GHz. Also shown are scatterplots of the same brightness temperature difference vs D 0 for 
various mass-density relationships at (b) 89 and (d) 150 GHz. 


7) At the CloudSat frequency of 94 GHz, radar back- 
scatter and attenuation are both very strongly sensi- 
tive to both D 0 and p for all but the smallest values of 
D 0g . Moreover, nearly solid ice particles with D 0 of 
2-3 mm — representative of large, dense graupel or 
sleet — exhibit two-way path attenuation as high as 
17 dB (kg m -2 ) -1 . These results would seem to pre- 
clude any snow water content retrieval from CloudSat 
alone without accurate independent knowledge of 
both p and D 0 , or constraints imposed by collocated 


passive microwave observations at 89 GHz and higher 
(e.g., from the scattering ratios from Fig. 7). 

4. Conclusions 

We have presented a simplified framework for assess- 
ing the qualitative sensitivities of computed microwave 
properties, satellite brightness temperatures, and radar 
reflectivities to assumptions about the physical properties 
of ice-phase hydrometeors. Properties considered included 
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a)S« Q /S, 


b ) S 89 /S 36.5 


10.0 


120 




p = 0.009 [g/cm 3 ] 
p = 0.018 [g/cm 3 ] 
p = 0.033 [g/cm 3 ] 

— 

p = 0.1 14 [g/cm 3 ] 

— 

p= 0.917 [g/cm 3 ] 


a =0.022, b =1.5 

o 

p p 

a =0.012, b =1.0 
p P 

+ 

a =0.015, b =1.18 

p P 

x 

a =0.018, b =0.8 

❖ 

p P 

a =0.015, b =0.6 

p p 


10.0 


const, p [g/cm ] 


C ) S 150 /S 89 


d ) S 150 /S 89 



0.18 


0.37 0.55 

const, p [g/cm 3 ] 


0.73 


0.92 


10.0 



Fig. 7. The two-frequency ratio of brightness temperature depression (relative to a 273-K background) for the finite ice layer described 
in Fig. 5. (a) The 89-GHz depression divided by the 36-GHz depression for constant-density PSDs. (b) As in (a), but for variable-density 
PSDs [with constant curves (solid lines)] for comparison, (c) The 150-GHz depression divided by the 89-GHz depression, (d) As in (c), but 
for variable-density PSDs. 


the shape and slope parameters of a gamma size distri- 
bution, particle shape and density, dielectric mixing 
formula, and the choice of complex index of refraction 
for ice. 

On the basis of our analysis, we find that all passive 
and active microwave observables of ice particles are 
extremely sensitive to the melted-equivalent mass me- 
dian diameter D 0 of the size distribution, all other pa- 
rameters being held constant. Similar large sensitivities 
are found for the ice fraction / ice whenever the geo- 
metric mass median diameter D 0g exceeds approxi- 
mately Vsth of the wavelength. Uncertainties in the 


correct specification of D 0 and/or / ice will likely domi- 
nate the errors associated with both forward model 
calculations and physical retrievals from satellite mi- 
crowave observations. 

By contrast, the distribution shape parameter /r has 
a surprisingly weak effect on any observable, amounting 
to no more than a 1-2-K variation in high-frequency 
brightness temperature for a 1 kg m -2 ice water path. 
For the radar backscatter at 13.4 and 35.6 GHz, the ef- 
fect of varying /jl from —1 to 4 was slightly more pro- 
nounced, yielding up to a 2.7-dB difference in the 
backscatter ratio for a distribution with D 0 = 1 mm. 



December 2012 


JOHNSON ET AL. 


2167 


Table 4. Results for selected radar frequencies. Baseline parameter values are indicated in boldface. Attenuation A is expressed as 

decibels per kilogram per meter squared of ice water path. 


13.4 GHz 35.6 GHz 94.0 GHz 



(dBZ) 

Z e (dBZ) 

2A (dB) 

Z e (dBZ) 

2A (dB) 

Ze (dBZ) 

2A (dB) 

/ ice (constant p) 

0.01 (0.0092) 

30.0 

25.1 

0.00 

12.4 

0.07 

-4.8 

0.58 

0.03 (0.028) 

30.2 

27.6 

0.01 

18.2 

0.11 

1.9 

1.13 

0.1 (0.092) 

30.5 

29.3 

0.01 

23.5 

0.19 

9.5 

2.35 

0.3 (0.28) 

30.8 

30.4 

0.01 

27.3 

0.30 

16.5 

5.04 

1.0 (0.92) 

30.2 

30.1 

0.01 

29.0 

0.37 

21.6 

10.42 

a p , b p [Eq. (2)] 

0.022, 1.50 

30.6 

29.6 

0.02 

25.4 

0.60 

18.4 

11.40 

0.012, 1.00 

30.5 

29.4 

0.01 

23.9 

0.21 

10.2 

2.65 

0.015, 1.18 

30.5 

29.4 

0.01 

24.0 

0.25 
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2.43 

Mixing rule 
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32.4 
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The choice of the Bruggeman versus Maxwell Garnett 
formulations for the effective dielectric constant of an 
ice-air mixture proved to be of at most very minor im- 
portance (<1.5 K and <0.3 dB), provided that one uses 
Maxwell Garnett with ice as the inclusion in an air matrix. 
Reversal of the roles of ice and air leads to a signifi- 
cant change in both microwave brightness tempera- 
ture (—10 K) and radar reflectivity (2 dB). 

For particle sizes and densities associated with strong 
scattering-based brightness temperature depressions at 
higher microwave frequencies (e.g., 89 and 150 GHz), 
switching from the older Warren (1984) index of refrac- 
tion to the newer Warren and Brandt (2008) produces 
a roughly 3%-4% increase in the magnitude of the 
scattering-induced depression when constant-density 
spheres are assumed. If one includes the Magono and 
Nakamura size-density relationship, then the difference 
between WB08 and W84 is as much as 14 K at 89 GHz 


and 27 K at 150 GHz. The choice has negligible effect 
on radar backscatter, however. 

The ratio of scattering-induced brightness temperature 
depressions for 36 and 89 GHz and for 89 and 150 GHz is 
found to be sensitive primarily to the melted-equivalent 
mass median diameter D 0 , with only very weak sensitivity 
to the density. These ratios, combined with the sensitivity 
of the GPM DPR to geometric mass median diameter 
D 0g , suggest that the DPR can directly complement GPM 
microwave imager (GMI) observations of the scattering 
depression to infer both particle size and density for 
a broad range of characteristic particle sizes, assuming 
one adopts one of the several available particle-density 
models. 

For the higher-frequency radars, the CloudSat radar 
at 94 GHz is extremely sensitive to both D 0 and p, and so 
retrievals of ice water content and/or precipitation ratio 
from this instrument alone would be extremely sensitive 
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Fig. 8. (left) The effective reflectivity factor Z e for W = lgm" 3 and (right) two-way path attenuation per 
1 kg m -2 of ice water path for 13.4-GHz (a) reflectivity and (b) attenuation, 35.6-GHz (c) reflectivity and 
(d) attenuation, and 94.0-GHz (e) reflectivity and (f) attenuation. 
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Fig. 9. Similar to Fig. 8, but including variable-density PSD -integrated reflectivities (points) and constant- 
density PSD-integrated reflectivities (solid lines). 
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to errors in either hydrometeor property. Moreover, 
two-way path attenuation per unit ice water path be- 
comes large (up to 18 dB) for dense ice particles with D 0 
larger than a few hundred micrometers. Although we 
did not present results for the highest frequencies of the 
Advanced Microwave Sounding Unit (183 GHz) or the 
GMI (166 GHz), our findings for 150 GHz are broadly 
applicable to those frequencies as well. 

In the interest of isolating the effects of variable ice 
particle properties, we did not consider the roles of 
cloud liquid water or gaseous absorption and emission 
on observed microwave brightness temperatures. It is 
understood that any additional absorber or attenuation 
in the atmospheric column will tend to reduce the overall 
sensitivity of microwave observables to parameter vari- 
ations we have described herein; we therefore present 
these results as “worst case” sensitivities. 

Although this study is theoretical in nature, it does lead 
to some recommendations to the broader remote sensing 
and validation communities who observe and measure 
ice-phase precipitation. The primary recommendation is 
to make direct measurements of the microwave extinc- 
tion (and scattering, if possible) by hydrometeors at 
several frequencies consistent with those employed on 
current passive and active microwave remote sensing 
platforms. This should also include coincident measure- 
ments of the physical properties of size, shape, fall speed, 
and so on such that our physical assumptions can be 
closely linked to the measured extinction properties. 
Having an instrumentation suite at various ground vali- 
dation sites that makes such measurements would be of 
significant benefit to reducing the above uncertainties. 
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